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Multigap superconductivity in sesquicarbides La 2 C3 and Y2C3 
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A complex structure of the superconducting order parameter in LnzCz {Ln = La, Y) is demon- 
strated by muon spin relaxation (/iSR) measurements in their mixed state. The muon depolarization 
rate [ov(T)] exhibits a characteristic temperature dependence that can be perfectly described by a 
phenomenological double-gap model for nodeless superconductivity. While the magnitude of two 
gaps is similar between La2C3 and Y2C3, a significant difference in the interband coupling between 
those two cases is clearly observed in the behavior of a v (T). 

PACS numbers: 74.70.Ad, 76.75.-K, 74.25.Jb 
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The revelation of high-temperature superconductivity 
in magnesium diboride (MgB2 , with critical temperature 
T c ~ 39 K) has stimulated renewed interest in other 
boride and carbide superconductors as an alternative 
path to novel superconductors with an even higher T c 
[y . Sesquicarbides (LTI2C3, Ln — La, Y) are among such 
compounds reported in early literatures; they exhibit su- 
perconductivity at relatively high critical temperatures 
(T c ~ 6-11 K) and their T c 's strongly depend on car- 
bon composition @, 0, 0]. Recently, we have found a 
new superconducting phase in Y2C3 that exhibits a much 
higher T c (~ 18 K) comparable with A-15 compounds 
Q . This discovery has attracted further attention to the 
relationship between structural details and superconduc- 
tivity in sesquicarbide systems. However, despite various 
attempts [6J, [7|, [8|, [9j , little is known so far about the de- 
tails of superconducting order parameters in La2C3 and 
Y2C3 from a microscopic viewpoint. 

A recent study on the temperature dependence of the 
nuclear spin-lattice relaxation rate in Y2C3 has suggested 
the occurrence of multiple superconducting gap with s- 
wave symmetry in a sample having T c — 15.7 K [icl |. 
While similar electronic structure would be expected for 
La2C3 0], a report on the specific heat measurement sug- 
gests single-gap superconductivity in a specimen with T c 
w 13.4 K [1JJ. In any case, the real nature of super- 
conductivity in Y2C3 and La2C3, including potential dif- 
ference between the two systems, still remains largely 
unclear. 

The muon spin rotation (/^SR) technique is a useful 
microscopic tool for probing quasiparticle (QP) density 
of states available for thermal/field-induced excitation in 
the mixed state of type II superconductors [H, Gj] • The 
muon depolarization rate in the mixed state is predomi- 
nantly determined by the magnetic penetration depth (A) 
that is controlled by superfluid density. Since the latter 
is reduced by the QP excitation, the effective value of A 
serves as a monitor of the QP excitation. In this letter, 
we present the result of /iSR measurements on polycrys- 
talline samples of La 2 C 3 (T c ~ 11 K) and Y 2 C 3 (T c ~ 15 
K), where a clear sign of double-gap superconductivity 
is observed in the temperature dependence of the muon 
depolarization rate. They also provide the first clear case 



for the double-gap model, where the magnitude of cou- 
pling between electronic bands responsible for supercon- 
ductivity is explicitly examined. Our result establishes a 
coherent description of multiple band/gap superconduc- 
tivity in this sesquicarbide system. 

For the Ln^Cs samples, starting materials were pre- 
pared by the arc melting method using a mixture of La/Y 
(99.9 %) and C (graphite, 99.99 %) with stoichiometric 
composition of sesquicarbide. The obtained Y-C alloys 
were placed into a BN cell in a dry box under an argon 
gas atmosphere, and polycrystalline Y2C3 was synthe- 
sized by elevating temperature to 1300 ~ 1400 °C for 30 
min under a high pressure of 5 GPa using a cubic-anvil- 
type equipment. For the polycrystalline La2C3, the La-C 
alloys obtained by the arc melting were pressed into pel- 
lets in a sealed Ta tube, and sintered at 1000 °C for 200 
h under a high vacuum condition of 3.0 x 10 -5 Torr, fol- 
lowed by a slow cooling process to ambient temperature 
at a rate of 5 °C/h. 

The powder x-ray diffraction patterns for both speci- 
mens could be indexed as a sesquicarbide phase with the 
space group of /43c?. In La2C3, nearly 10% of LaC2 was 
observed as a minor phase besides that of the sesquicar- 
bide, while Y2C3 was found to be in a single phase. LaC2 
behaves as a normal metal above 2 K and only causes 
a background in the /iSR signal in the superconducting 
phase. The lattice constants of La2C3 and Y2C3 were 
determined to be approximately a — 8.808(5) A and 
8.238(5) A, respectively, which are in good agreement 
with those reported previously [III EH EE EB| ■ Unfortu- 
nately, the precise stoichiometry of carbon has not been 
determined. Therefore, the chemical composition in this 
paper refers only to a nominal value. Heat capacity and 
ac and dc magnetic susceptibilities were measured using 
MPMSR2 and PPMS (Quantum Design Co., Ltd.). 

Figure [1] shows the temperature dependence of dc sus- 
ceptibility under zero- field cooling in La2C3 and Y2C3. 
In both cases, clear diamagnetic signals are observed 
below T c (~ 11.2 K and 15.2 K, respectively). As 
shown in the inset, H C 2(T) exhibits almost linear depen- 
dence on temperature, which differs significantly from 
the Werthamer-Helfand-Hohenberg (WHH) relation [lTj ■ 
The enhancement of H 2(T) from the WHH prediction 



4 8 12 16 

T(K) 

FIG. 1: (Color online) Temperature dependence of dc mag- 
netic susceptibility at 10 Oe in La2C3 and Y2C3 normalized 
by the value at 2 K. Inset shows magnetic field (H) vs. tem- 
perature (T) phase diagram. The triangle, circle, and square 
symbols indicate data determined by heat capacity and dc and 
ac magnetic susceptibility measurements, respectively. The 
open and closed symbols show data obtained from T and H- 
scan, respectively. The solid and dashed curves correspond to 
the WHH relation and the GL model, respectively. 



is due to the strong electron-phonon coupling rather 
than the anisotropic Fermi surface or localization ef- 
fect [ll[. We can extract H c2 (0) without much uncer- 
tainty using the Ginzburg-Landau (GL) theory, H C 2(T) 
= H c2 (0)(l - (T/T c ) 2 )/(1 + (T/T c ) 2 ), where H c2 {0) = 
$o/27t£gl i ^0 is the flux quantum, and £gl is the GL- 
coherence length . The best fit using the above equa- 
tion yields H c2 (0) = 167(3) and 256(7) kOe for La 2 C 3 
and Y2C3, respectively. 

Conventional /iSR experiment was performed on the 
M15 beamline of TRIUMF, Canada. The polycrystalline 
samples were loaded on a sample holder (a scintillator 
serving as a muon veto counter, with a sample dimension 
of 7 x 7 mm 2 ) and placed into a He gas- flow cryostat, to 
which a 100% spin-polarized muon beam with a momen- 
tum of 29 MeV/c was irradiated to collect 1.5 x 10 7 decay 
positron events for each spectrum (taking about 1.5 h). 
Each measurement was performed under a field-cooling 
process to minimize the effect of flux pinning, and field 
fluctuation was kept within 10 -4 of the applied field. 

Since we can reasonably assume that muons stop ran- 
domly on the length scale of the flux- line lattice (FLL), 
the muon spin precession signal, P(t), provides the ran- 
dom sampling of the internal field distribution B(r), 

/oo 
n(B) cos^Bt + </>)dB, (1) 
-OO 

n(B) = (S(B-B(r))) r (2) 

where 7 M is the muon gyromagnetic ratio (= 27r x 13.553 
MHz/kOe), n(B) is the spectral density for the inter- 
nal field defined as a spatial average (() r ) of the delta 
function, and <f> is the initial phase of rotation. These 
equations indicate that the real amplitude of the Fourier 
transformed muon spin precession signal corresponds to 
n(-B) (except corrections for additional relaxation due 
to other origins, see below). In the case of relatively 



FIG. 2: (Color online) Time revolution of muon-positron de- 
cay asymmetry in (a) La2C3 and (b) Y2C3 at 2 K under a 
transverse field of 2.5 kOe and 5.0 kOe, respectively, displayed 
in a rotating-reference-frame frequency of (a) 33 MHz and (b) 
66.8 MHz. The respective insets show the fast Fourier trans- 
form (FFT) at 2 K (solid lines) and above T c (dashed lines). 



large magnetic penetration depth (A > 3000 A), n(B) 
can be well-approximated by a simple Gaussian field pro- 
file, yielding P(t) ~ exp(— a 2 t 2 /2) cos^^Bt + </>), where 

a — 7 M \/((-B — B( r )) 2 ) fx ^ 2 an d B ~ H is the mean 
field. Here, it must be stressed that A is an effective mag- 
netic penetration depth susceptible to the quasiparticlc 
excitation. 

Figure [2] shows the time-dependent muon-positron de- 
cay asymmetry at 2 K in La2C3 and Y2C3 with their 
fast Fourier transform (FFT) displayed in the inset. The 
FFT spectral linewidth in the normal state (T > T c ) is 
determined by the small random local fields from nuclear 
moments and a limited /iSR time window (~ 8 /xs), while 
that in the superconducting state is further broadened by 
the formation of FLL and associated inhomogeneous lo- 
cal field distribution [-B(r)]. The solid curves in the main 
panels are the best fits of the data in the time domain, 
assuming two components of the Gaussian damping, 

2 / CT 2£2 \ 

AP{t) = ^A.expf — *— J cos( 7pJ B 4 < + <fc) (3) 

i=l ^ ' 

where the i-th component refers to the contribution from 
superconducting (i — 1) and normal (i — 2) phases, A4 is 
the partial asymmetry (Y]. Aj — A), ai is the relaxation 
rate, and 7 M -B» is the central frequency for the respec- 
tive components. The model yields good fits to data, 
as indicated by the reasonably small values of reduced 
chi-square: % 2 /Nf is mostly less than 1.7 for La2C3 and 
1.3 for Y2C3, with Nf being the number of degrees of 
freedom. Considering that a 2 represents the relaxation 
due to the nuclear magnetic moment (i.e., a 2 — a n ), the 
net relaxation rate in the superconducting state is ex- 
pressed as a\ = a 2 + a 2 , where the second term comes 
from n(B) in the FLL state and it is proportional to the 
superfluid density From the fitting analysis, the su- 
perconducting volume fractions [= Ai/(Ai +^4 2 )] at 2 K 
are estimated to be w 0.91 and 0.98 in La 2 C3 and Y 2 C3, 
respectively, where the former value is in good agreement 
with the fractional yield estimated by the X-ray analy- 
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FIG. 3: (Color online) Temperature dependence of muon spin 
relaxation rate for (a) La2C3 at 2.5 kOe and (b) Y2C3 at 5.0 
kOe. Error bars (not shown) are smaller than the symbol 
size. Solid and dashed curves indicate the result of fitting 
analysis using the double-gap model described in the text. 
Insets show the relaxation rate in the superconducting state 
(<Ti) as a function of magnetic field for La2C3 (a) and the 
order parameters [A(T)//sbTc] f° r the respective cases. 



sis (with the rest corresponding to LaC2 in the normal 
state). 

Figure [3] shows the temperature dependence of cr v for 
La2C3 at 2.5 kOe and Y2C3 at 5.0 kOe, where the fields 
are chosen to be reasonably away from the lower critical 
field (the field dependence at the relevant field range is 
shown for L^Cs). It is interesting to note that the data 
of La2C3 exhibit a shoulder-like structure near 4 K, in- 
dicating clear deviation from the behavior predicted for 
single-gap BCS superconductors. The effect of flux pin- 
ning as a possible origin of such structure is ruled out 
by the fact that o\ at 2 K is almost independent of the 
applied field below 30 kOe [see the inset of Fig. [3Ja)]. 
According to the theories that consider multiband super- 
conductivity [13, HH, H3] , such an inflection is expected to 
occur in the systems consisting of two superconducting 
bands that are weakly coupled. The behavior of er v below 
^4 K is then attributed to the band with a smaller gap 
energy. To our knowledge, this is the first unambiguous 
example of double-gap superconductivity with extremely 
week interband coupling. 

Compared with La2C3, no strong anomaly is observed 
in the case of Y2C3. However, while the superfluid den- 
sity (and hence er v ) is predicted to be virtually indepen- 
dent of temperature for T/T c < 0.4 in the single-gap BCS 
model, a clear variation of <r v with temperature is ob- 
served below 6 K. This is another sign that the order 
parameter of Y2C3 has an anisotropic structure. Con- 
sidering the present result in La2C3 and previous NMR 
studies [l(| , we can attribute the T-dependence of a v also 
to the double-gap superconductivity. This is further sup- 
ported by a recent first-principles calculation [23| , where 
it is suggested that the Fermi surface of Y2C3 consists of 
three hole bands (0.12, 0.15, and 0.88/eV unit cell spin) 
and one electron band (2.73/eV unit cell spin) that arise 



TABLE I: Superconducting properties of La2C3 and Y2C3 de- 
termined from the present experiment, where those obtained 
from the double-gap analysis correspond to the solid curves 
in Fig. |3 









Transverse field (kOe) 


2.5 


5.0 


T c (K) 


10.9(1) 


14.7(2) 


o-v(O) (fis' 1 ) 


0.71(3) 


0.48(2) 


A(0) (A) 


3800(100) 


4600(100) 


10 


0.38(2) 


0.86(2) 


Ai(0) (meV) 


2.7(1) 


3.1(1) 


A 2 (0) (meV) 


0.6(1) 


0.7(3) 


2A 1 /fc B T c 


5.6(3) 


4.9(3) 


2A 2 /fc B T c 


1.3(3) 


1.1(5) 



mainly from the hybridized orbitals of Y d- and C-C an- 
tibonding 7r*-states. Such large differences in the density 
of states and Fermi velocities between hole and electron 
bands might lead to the opening of two superconducting 
gaps in the different parts of the Fermi surface. 

The origin of difference in the temperature dependence 
of cr v between La2C3 and Y2C3 is understood by consid- 
ering the difference in the interband coupling strength be- 
tween these two compounds. For quantitative discussion, 
the data in Fig. [3] are analyzed using a phenomenological 
double-gap model with s-wave symmetry [2~H l25j , 

<r(T) = tr(Q) - w5a(Ai,T) - (1 - w)-Sa(A 2 ,T), 

6a(A,T) = ^J™f(e,T)-[l-f(e,T)]de, 

f(e,T) = (l + eV £2+A(T)2/fcBT )~\ 

where Aj (i = 1 and 2) is the energy gap at T = 0, 
w is the relative weight for i = 1, fes is the Boltzmann 
constant, /(e, T) is the Fermi distribution function, and 
A(T) is the standard BCS gap energy. The solid curves 
in Fig. [3] are the best fit result obtained by using the 
above double-gap model with the parameters listed in 
Table HI For La2C3, a simplified model (dashed curve) 
assuming two independent superconducting bands was 
also tested against the data, which turned out to exhibit 
a slightly better agreement than that described by us- 
ing the above model [yielding 2/S.i/k^,T c = 4.5(3) and 
2A 2 /fcBT c = 1.3(3)]. This might suggest that the above 
model may not necessarily be a good approximation for 
the case of weak interband coupling. 

The superconducting parameters deduced from the 
present experiment are summarized in Table |TJ Here, we 
calculated the magnetic penetration depth A(0) using the 
formulaCT v (0) [fis" 1 } = 4,83xlOHl-H/H c2 )x [1+3.9(1- 
H/H c2 ) 2 ] 1/2 /\ 2 {Q) [nm] [H H|. The gap parameter 
2A/fce7c of Y2C3 is in reasonable agreement with that 
deduced by NMR (i.e., 2A/fc B T c = 5 and 2) 10], again 
supporting the present double-gap scenario. We also find 
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that 2A//cbT' c for the two respective bands of La2C3 are 
comparable with those of Y2C3. Thus, it appears that 
the superconductivity of La2C3 and Y2C3 share the com- 
mon features of strong electron-phonon coupling and s- 
wave symmetry, which is in line with the previous heat 
capacity results fill ]. 

Provided that there is a significant difference in the 
interband coupling between La2C3 and Y2C3, the ob- 
served difference in the relative weight (w) between two 
gaps might also be connected with the interband cou- 
pling. Furthermore, considering that the double-gap fea- 
tures tend to be suppressed by the localization (scatter- 
ing) effect, one might suspect that such a difference in w 
may arise from that in the quality of specimen. In this 
regard, we have to note that the Y2C3 samples were ob- 
tained only in a polycrystalline form using high pressure 
synthesis and that their short annealing time might have 
resulted in a quality less than that of La2C3. At this 
stage, we presume it unlikely that the present result has 
been strongly affected by the localization effect, consider- 
ing that the electronic mean free path measured using the 
microwave cavity perturbation technique is much longer 
than £gl for the sample prepared under the same condi- 
tion [23. However, it would be certainly helpful to study 
the influence of sample quality (and chemical stoichiom- 
etry as well) in the future to elucidate the details of the 
localization effect on the double-gap behavior. 

Finally, let us point out the noncentrosymmetric ef- 
fect in superconductivity. In the case of a sesquicar- 
bides system with the /43d group symmetry, an asym- 
metric spin-orbit interaction can be approximated by the 
Dresselhaus-type interaction. When the order of magni- 



tude of a superconducting gap is of comparable to that 
of the spin-orbit band splitting, the original isotropic 
gap structure is modulated by a magnetic field to have 
a point-node, because anisotropic Pauli depairing effect 
can occur in the specific part of the momentum space 
[2^ |. This may lead to unusual field-induced quasipar- 
ticle excitation, and a detailed /iSR study on the field 
dependence of magnetic penetration depth is currently 
in progress to examine the proposed scenario. 

In summary, we performed /iSR experiment on L712C3 
(Ln = La, Y) to clarify the structure of superconducting 
order parameter through the temperature dependence of 
quasiparticle excitation reflected in the muon depolar- 
ization rate, cr v (T) in the mixed state. We showed that 
<J V (T) exhibits a characteristic of double-gap in the su- 
perconducting order parameter, with a marked variation 
in the temperature dependence between La and Y com- 
pounds that is attributed to the difference in the inter- 
band coupling. The gap parameters for two respective 
bands were deduced using the phenomenological double- 
gap model and were found to be comparable between La 
and Y compounds, which is consistent with the occur- 
rence of a strong-coupling superconductivity with s-wave 
symmetry in both the systems. 
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